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Ultracold Fermi atoms allow the realization of the crossover from 
Bardeen-Cooper-Schrieffer (BCS) superconductivity to Bose- Ein- 
stein condensation (BEC) 1 3 , by varying with continuity the attrac- 
tion between fermions of different species 4 . In this context, radio 
frequency (RF) spectroscopy 5 7 provides a microscopic probe to 
infer the nature of fermionic pairing. In the strongly-interacting 
regime, this pairing affects a wide temperature range comprising 
the critical temperature T c , in analogy to the pseudogap physics for 
high-temperature superconductors. By including what are known 
in superconductors as "paraconductivity" fluctuations, here we cal- 
culate the RF spectra above T c for fermions with balanced popula- 
tions and compare them with available experimental data, thus re- 
vealing that the role of these fluctuations is considerably enhanced 
with respect to superconductors. In addition, we extract from the 
spectra an energy scale associated with pairing and relate it to a 
universal quantity recently introduced for Fermi gases 8 . 
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It is widely appreciated at present that ultracold Fermi atoms represent 
ideal systems for deepening our understanding of condensed-matter systems, 
especially as far as the many-body quantum physics is concerned. This is 
because the inter-particle interaction and species populations can be fine- 
controlled and tuned almost at will, in such a way that the system Hamil- 
tonian is precisely known over a wide parameter range. Since experimental 
techniques are able to provide diverse and accurate data, major theoretical 
efforts can be justified to attempt detailed comparison with experiments, 
particularly for those aspects that have remained unsettled in condensed 
matter. 

In this respect, RF experiments 5-7 (especially in their tomographic ver- 
sion) provide fertile ground for exploring complementary aspects related to 
excitation gaps and pairing fluctuations, both below and above the critical 
temperature T c for the transition to the superfluid phase in these neutral 
systems. A crucial aspect here is the occurrence of several hyperfine lev- 
els, which are split by a magnetic field and exhibit mutual Fano-Feshbach 
resonances that strongly amplify the two-body interaction when the field is 
varied 4 . The hyperfine levels act as components of a (hyper) spin and a spin 
flip is generated by the RF transition, so that a dynamical (two-particle) 
spin-correlation function is effectively measured by RF spectra 9 . 

Below T c , theoretical interpretation of these spectra has revealed a com- 
petition between pairing-gap effects in the initial state of the transition and 
final-state effects, which tend to push the oscillator strength toward opposite 
sides of the spectrum 10 . In this case, the presence of a well-developed pairing 
gap at low temperature makes a BCS-RPA approach 11 for response functions 
sufficiently accurate to reproduce the experimental features 10 . 

When the temperature increases across T c the pairing gap gives the way 
to pairing fluctuations, and precursor effects of pairing appear in the normal 
phase above T c . These effects manifest themselves at the simplest level in 
the single-particle spectral function with a depression of the spectral weight 
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(pseudogap) at low frequencies for small enough wave vectors (the effect being 
amplified when approaching the BEC limit of the BCS-BEC crossover 12 ). In 
high-temperature cuprate superconductors (where pairing is competing with 
other kinds of ordering), photoemission (ARPES) experiments have exten- 
sively been used to unravel the nature of the pseudogap in single-particle 
excitations 13 . Recently, attempts to disentangle single-particle excitations 
from RF (two-particle) spectra have also been made for ultracold Fermi 
gases 14 , by resolving the wave vector of the photo-excited fermions. 

Pairing fluctuations in the normal phase above T c are known to affect the 
two-particle response of superconductors. In particular, for superconducting 
thin films the possibility of fluctuational creation of Cooper pairs has been 
shown to enhance their conductivity, an effect known as "paraconductivity" 
from the work of Aslamazov and Larkin 15 (AL). Another source of inter- 
action of fermions with fluctuating Cooper pairs stems from the additional 
contribution to the conductivity found by Maki 16 and Thompson 17 (MT). 
The AL and MT contributions have been detected experimentally to pro- 
duce small changes to the normal conductivity, which can be amplified by 
reducing the dimensionality 18 . While the AL contribution to the conductiv- 
ity can be obtained theoretically by a time- dependent generalization of the 
Ginzburg-Landau equations 19 , accounting for the MT contribution requires 
one to introduce finite-temperature Feynman diagrams for the current re- 
sponse function, a technique which yields also a third relevant contribution 
known as density-of-states (DOS) renormalization 20 . 

In this paper, we calculate the corresponding fluctuation contributions 
to the dynamical spin-correlation function associated with RF spectra, by 
using finite-temperature Feynman diagrams in the normal phase. Our calcu- 
lation reveals that, even in three dimensions, the role of pairing fluctuations 
can be drastically enhanced for the RF spin-correlation function of ultracold 
Fermi atoms with respect to the current-correlation function of superconduc- 
tors. The results of our calculation compare favorably with the available RF 
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experimental data for Fermi gases with balanced populations under various 
coupling conditions. 

The differences between the current- and RF spin-correlation functions 
originate from the different nature of the coupling to the external probe, the 
involvement of a third spin component and the wider coupling range that 
can be explored via Fano-Feshbach resonances with Fermi atoms. Apart 
from these differences (which yet make the MT contribution vanish identi- 
cally for the RF spin-correlation function), the topological structure of the 
Feynman diagrams for the AL and DOS contributions remains the same for 
both calculations (see the Methods section). 

For the RF spectra of interest here, two spin components (a and f3) are 
initially equally populated while a third (7) is empty. For a transition which 
flips (3 into 7, the relevant scattering lengths are then a a p for the initial (i) 
and a ai for the final (/) state. In our calculation, the DOS contribution 
contains only pairing fluctuations related to a al 3 (initial-state effects), while 
the AL contribution contains in addition pairing fluctuations related to a ai 
(final-state effects) - see the Methods section. For this reason, we refer to 
a a/ 3 as di and to a ai as a/. Experiments have been able to explore differ- 
ent combinations of these spin components, thereby enlarging the accessible 
coupling ranges. Figure 1 reports the corresponding locations of the possible 
experimental coupling values for 6 Li in the ((/cf^) -1 , (/c^a/) -1 ) plane, where 
&f is the Fermi wave vector. 

Comparisons between the results of our numerical calculations and the 
experimental spectra for 6 Li are reported in Fig. 2, resulting in quite good 
agreement between theory and experiments. Recall that, in the absence of 
pairing fluctuations, a single delta-spike at zero frequency would absorb the 
whole spectra weight (whose value is dictated by a sum rule for the spec- 
tral area 10 and is not changed by further inclusion of a« and a/). In all 
panels, the left peak is associated with the presence of a localized (bound) 
final state, which typically absorbs a spectral weight comparable to the right 
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peak. The latter originates from continuum transitions and comprises the 
high-frequency tail. For these reasons, the final-state effects associated with 
the AL contribution play an essential role in all cases (see also Supplementary 
Fig. 1). Nonetheless, when the bound and continuum peaks are sufficiently 
apart, the continuum part of the spectra appears to be only slightly modified 
by the inclusion of final-state effects provided the peak heights (obtained by 
DOS only and by DOS plus AL) are brought to coincide. These considera- 
tions justify a posteriori the procedure adopted in ref. 7 to fit the continuum 
part of the experimental spectra disregarding final-state effects. 

Experimentally, the motivation for exploring RF spectroscopy with ul- 
tracold Fermi atoms 5 stemmed from the expectation of extracting the value 
of a pairing gap directly from the (continuum part of the) spectra. This 
would be possible within a BCS description al low temperature, since in this 
case the continuum edge occurs at A 2 / (2E F ) where A is the BCS gap pa- 
rameter and E F = k F /(2m) the Fermi energy (m being the fermion mass). 
Such a simplified picture, however, does not apply to the unitarity regime 
about {kFO,i)~ l = owing to strong-coupling effects in the initial and/or final 
states of the transition. Accordingly, no procedure has been suggested so far 
for extracting from the RF spectra an energy scale associated with pairing 
interaction in systems with balanced populations 23 . 

The present analysis allows us to identify such an energy scale when the 
bound and continuum peaks of an RF spectrum are sufficiently far apart 
(corresponding to the green area of Fig. 1). In this case, there exists a win- 
dow on the high-frequency side of the continuum peak where the RF signal 
behaves like B(lo/E f )- 3 / 2 with B = A (3/2 5/2 ) (A^/Ep) 2 . Here, A stands 
for the area of the continuum peak once the total area of the RF spectrum 
is taken to be unity. Typically, this frequency window extends between a 
few times Ep and (ma 2 ) -1 , past which the spectral tail is dominated by 
final-state effects and decays like to~ 5 ^ 2 . The above expression is written in 
analogy with the high-frequency behavior of the RF signal within the BCS 
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approximation, in which A appears in place of A^. The energy scale Aoo 
can be further related to the asymptotic behavior C/k 4 of the wave- vector 
distribution function n(k), where C is the "contact intensity" introduced in 
ref. 8 which enters several quantities of a Fermi gas in a universal way 24 ' 25 . 
From our theory we identify C = (mA^) 2 (see the Methods section). 

In Fig. 3a the theoretical RF spectra show a plateau whose height is 
(3/2 5 / 2 ) (Aoo/Ep) 2 according to the above argument. This height is com- 
pared with the horizontal line, obtained by the value of A^/Ep resulting 
from our independent calculation of the contact intensity C through n(k). 
In Fig. 3b we apply the same procedure to the experimental RF spectra. Al- 
though the data are somewhat scattered in the frequency window of interest, 
their tendency to approach a plateau results from the figure and allows us to 
obtain a numerical estimate for A^. This procedure could be implemented 
in future experiments to map the value of A^ in extended temperature and 
coupling ranges. 
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METHODS 



RF SPECTRUM AND THE SPIN CORRELATION FUNCTION 

The RF spectrum is obtained theoretically within linear-response theory as: 

8(I(u th )) = -2g 2 jdvdv' Im{n*(r,r> th )} 

where g is the coupling constant of the atomic {(5 — > 7) transition, Hf 7 (r, r'; c^th) 
is the Fourier transform of the retarded correlation function Hf 7 (r, r'; t—t') = 
—iO{t — t / )([5 / 3 7 (r, t), -Bg 7 (r', t')]) taken at the shifted frequency cu t h = ujrf + 
Hp — /x 7 with respect to the RF frequency ujrf (we set h — 1 throughout). 
Here, Bp (r) = ^(r^^r) is the transition operator in terms of the field 
operators ^g/ 7 at the spatial point r, /ip and /x 7 are the chemical potentials 
for the initial (equally populated) and final (empty) levels, in the order, and 
(• • ■) stands for a thermal average. 

Following a standard procedure 26 , the retarded correlation function is 
conveniently calculated via its Matsubara counterpart 

rWr,r';0,) = J* d r e^(T T [^(r', 0)^(r, r + )^ 7 (r, r)^(r', + )]) , 

for this expression admits a representation in terms of Feynman diagrams. 
Here, T r is the time-ordering operator for imaginary time r and Q u = li^vT 
[y integer) is a bosonic Matsubara frequency at temperature T (we set ks — 1 
throughout). Above T c , the relevant Feynman diagrams for the two-particle 
response are the DOS and AL contributions depicted in Fig. 4, which result 
once the single-particle self-energy with pairing fluctuations is adopted. An- 
alytic continuation from ifl u to cj t h + if] (v being a positive infinitesimal) is 
eventually required. To compare with the experimental spectra, we have to 
convert from u t h to u = Urf — ^7/3 where u) 1 p is the Bohr frequency of the 
atomic transition. Besides being appropriate to the BCS ((aapkp) -1 ~ —1) 
and unitarity ((cta/jA;^)™ 1 m 0) regimes, the DOS and AL contributions are 
also able to recover the two-body (molecular) calculation of ref. 27 when 
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the BEC regime (1 ~ (a a /3^F) x ) is approached. (In the text, we have set 
NUMERICAL PROCEDURES 

Due to the involved structure of the AL contribution (see Fig. 4b), for the ex- 
ternal Matsubara frequency Q u analytic continuation from the points ifl u on 
the imaginary axis to (just above) the real frequency axis cannot be achieved 
by the elementary methods that proved sufficient 10 for the BCS-RPA ap- 
proach below T c . The delicate analytic continuation at finite temperature 
which is required for the AL diagram has been performed in the literature 
with the resort to approximations 20 which do not apply to the present case, 
when experimental RF spectra with nontrivial frequency structures need be 
accounted for in a variety of situations. 

To this end, we have performed the analytic continuation of the AL con- 
tribution resorting to Pade approximants 28 , whereby the dependence on the 
external frequency Q u is approximated by the ratio of two polynomials (the 
denominator having one additional power with respect to the power M of the 
numerator). The 2M unknown coefficients have been determined by sam- 
pling over about 450 Matsubara frequencies (not necessarily equally spaced). 
The reliability of the numerical results has been tested against the DOS con- 
tribution above T c as well as the BCS-RPA calculation below T c , which can 
both be obtained by direct analytic continuation without resorting to Pade 
approximants. Additional independent tests on the numerical procedure are 
the values of the sum rules for the spectral area 10 and its first moment 9 , 
which are reproduced within 1% and 10%, in the order. 

This procedure of analytic continuation was avoided in ref. 29 while 
addressing the calculation of RF spectra above T c , by effectively extending 
the BCS-RPA calculation to temperatures between the actual T c and the 
corresponding temperature obtained at the mean-field level (often referred 
to as the temperature T* at which pairing exhausts its effects). 
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THE ENERGY SCALE Aoo OF THE PAIRING INTERACTION 

The most direct way to introduce the quantity is through the coefficient 
(3/2 5 / 2 ) A 2 ^ of the u~ 3 ^ 2 behavior of the DOS and AL contributions for large 
ou, which identifies 



where q is a wave vector and T a p is the pairing propagator of Fig. 4c. This 
expression is conveniently evaluated by using the spectral representation 30 
of F a p. In the BCS regime A^ coincides with the absolute value 2n\ai\n/m 
of the mean-field shift for temperatures smaller than (ma 2 ) -1 , while in the 
BEC regime the relation A 2 ^ = Aim/ (m 2 <ii) holds in analogy with the BCS 
mean-field result. Close to T c at unitarity, the numerical value 0.75E F of 
Aqo is comparable with the value 0.8E F of the pseudogap extracted from the 
single-particle spectral function 12 . The temperature dependence of A^ is 
rather weak for all coupling regimes, extending up to several times the Fermi 
temperature T F . 

The energy scale A^ is also related to the asymptotic behavior of the 
wave-vector distribution function n(k) as obtained from the fermionic single- 
particle propagator (3- (3 dressed by the self-energy of Fig. 4c. We get 
n(k) = (m Aoo) 2 /k 4 where Aj^ is defined as above. From this expression we 
identify (m A^) 2 with the contact intensity C of ref. 8. 




9 



References 

1. Eagles, D. M. Possible pairing without superconductivity at low carrier 
concentrations in bulk and thin-film superconducting semiconductors. Phys. 
Rev. 186, 456-463 (1969). 

2. Leggett, A. J. in Modern Trends in the Theory of Condensed Matter (Proc. 
16th Karpacz Winter School Theor. Phys.) (eds Pekalski, A. & Przystawa, 
J.) 13-27 (Springer, 1980). 

3. Nozieres, P. & Schmitt-Rink, S. Bose condensation in an attractive fermion 
gas: from weak to strong coupling superconductivity. J. Low Temp. Phys. 
59, 195-211 (1985). 

4. See, e.g., Ketterle, W. & Zwierlein, M. W. in Ultra- Cold Fermi Gases 
(Proc. Internat. School Phys. Enrico Fermi, Course 164) (eds Inguscio, M., 
Ketterle, W. & Salomon, C.) 95-287 (IOS Press, 2008). 

5. Chin, C. et al. Observation of the pairing gap in a strongly interacting 
Fermi gas. Science 305, 1128-1130 (2004). 

6. Shin, Y., Schunck, C. H., Schirotzek, A. & Ketterle, W. Tomographic 
rf spectroscopy of a trapped Fermi gas at unitarity. Phys. Rev. Lett. 99, 
090403 (2007). 

7. Schunck, C. H., Shin, Y., Schirotzek, A. & Ketterle, W. Determination 
of the fermion pair size in a resonantly interacting superfluid. Nature 454, 
739-743 (2008). 

8. Tan, S. Large momentum part of a strongly correlated Fermi gas. Ann. 
Phys. doi:10.1016/j.aop.2008.03.005 (2008). 

9. Yu, Z. & Baym, G. Spin-correlation functions in ultracold paired atomic- 
fermion systems: Sum rules, self-consistent approximations, and mean fields. 
Phys. Rev. A 73, 063601 (2006). 

10. Perali, A., Pieri, P. & Strinati, G. C. Competition between final-state and 
pairing-gap effects in the radio-frequency spectra of ultracold Fermi atoms. 
Phys. Rev. Lett. 100, 010402 (2008). 

11. Schrieffer, J. R. Theory of Superconductivity Ch. 8 (W. A. Benjamin, 

10 



New York, 1964). 

12. Perali, A., Pieri, A., Strinati, G. C. & Castellani, C. Pseudogap and spec- 
tral function from superconducting fluctuations to the bosonic limit. Phys. 
Rev. B 66, 024510 (2002). 

13. Damascelli, A. Probing the electronic structure of complex systems by 
ARPES. Phys. Scr. T109, 61-74 (2004). 

14. Stewart, J. T., Gaebler, J. P. & Jin, D. S. Using photoemission spec- 
troscopy to probe a strongly interacting Fermi gas. Nature 454, 744-747 
(2008). 

15. Aslamazov, L. G. & Larkin, A. I. Effect of fluctuations on the properties 
of a superconductor above the critical temperature. Soviet Solid State 10, 
875-880 (1968). 

16. Maki, K. Critical fluctuation of the order parameter in a superconductor. 
Prog. Theor. Phys. (Kyoto) 40, 193-200 (1968). 

17. Thompson, R. S. Microwave, flux flow, and fluctuation resistance of dirty 
type-II superconductors. Phys. Rev. B 1, 327-333 (1970). 

18. Skocpol, W. J. & Tinkham, M. Fluctuations near superconducting phase 
transitions. Rep. Prog. Phys. 38, 1049-1097 (1975). 

19. Abrikosov, A. A. Fundamentals of the Theory of Metals Ch. 19 (North- 
Holland, Amsterdam, 1988). 

20. Varlamov, A. & Larkin, A. Theory of Fluctuations in Superconductors 
Ch. 7 (Oxford Scholarship Online, Oxford, 2007). 

21. Bartenstein, M. et al. Precise determination of 6 Li cold collision parame- 
ters by radio-frequency spectroscopy on weakly bound molecules. Phys. Rev. 
Lett. 94, 103201 (2005). 

22. Simonucci, S., Pieri, P. & Strinati, G. C. Broad vs. narrow Fano-Feshbach 
resonances in the BCS-BEC crossover with trapped Fermi atoms. Europhys. 
Lett. 69, 713-718 (2005). 

23. Schirotzek, A., Shin, Y., Schunck, C. H. & Ketterle, W. Determination 
of the superfluid gap in atomic Fermi gases by quasiparticle spectroscopy. 



11 



Preprint at Qhttp://arxiv.org/abs/0808.0026v2D (2008). 

24. Braaten, E. & Platter L. Exact relations for a strongly interacting Fermi 
gas from the operator product expansion. Phys. Rev. Lett. 100, 205301 
(2008). 

25. Werner, F., Tarruell L. & Castin, Y. Number of closed-channel molecules 

in the BEC-BCS crossover. Preprint at ( |http://arxiv.org/abs/0807.0078vl[ )(2008). 

26. Fetter, A. L. & Walecka, J. D. Quantum Theory of Many-Particle Sys- 
tems Ch. 9 (McGraw-Hill, New York, 1971). 

27. Chin, C. & Julienne P. S. Radio-frequency transitions on weakly bound 
ultracold molecules. Phys. Rev. A 71, 012713 (2005). 

28. Vidberg, H. J. & Serene, J. W. Solving the Eliashberg equations by 
means of N— point Pade approximants. J. of Low. Temp. Phys. 29, 179- 
192 (1977). 

29. He, Y., Chien, C. C, Chen, Q. & Levin K. Temperature and final state 
effects in radio frequency spectroscopy experiments on atomic Fermi gases. 
Preprint at ( |http: / /arxiv.org/abs/08041429vl[ ) (2008) . 

30. Pieri, P., Pisani, L. & Strinati, G. C. BCS-BEC crossover at finite tem- 
perature in the broken-symmetry phase. Phys. Rev. B 70, 094508 (2004). 

Acknowledgements 

This work was partially supported by the Italian MUR under contract PRIN- 
2007 "Ultracold Atoms and Novel Quantum Phases" . This paper is dedicated 
to the memory of Professor Franco Bassani. 



12 



-2 

-1.5 -1 -0.5 0.5 1 1.5 

(k F aj)" 1 

Figure 1: Allowed coupling values for 6 Li in the ((k F ai) _1 , (k F a f ) ~ 1 ) 
plane. RF experiments realized with ultracold 6 Li atoms correspond to the 
three following combinations of scattering lengths: (I) for (a-i = ai2,«/ = 
013), (II) for (a« = 013, a/ = au), and (III) for (a« = ai 3 , a/ = 023). Using the 
values of these scattering lengths vs. the magnetic field given in ref. 21 and 
the experimentally accessible values of the Fermi wave vector kp = (37r 2 n) 1 / 3 
where n is the (total) density at the trap center (typically, kp ~ 2.7 x 10~ 4 ciq 1 
where a is Bohr radius), we construct the three stripes labeled (I), (II), and 
(III) inside which RF spectra can be collected, whose widths correspond to 
an estimated 40% variation on the values of kF- The symbols correspond to 
the three different sets of tomographic RF spectra reported in ref. 7. The 
green area indicates the region 2 ~ |(/c^aj) _1 — (fc^a^) _1 | where it is possible 
to extract the quantity A ro (see below) from the high-frequency tail of the 
RF spectra. 
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Figure 2: Comparison between theoretical and experimental RF spectra. 

The experimental spectra from Figs. 4a-d of ref. 7 (black squares) are compared 
with our theoretical calculations, obtained with the inclusion of pairing fluctuations 
due to di (DOS contribution - red dotted curves) and with the additional inclusion 
of pairing fluctuations due to aj (DOS plus AL contributions - blue full curves). 
The values aj, aj, and kp are taken from the experiments. The experimental 
spectra (which are tomographic for panels a-c and trap-averaged for panel d) span 
approximately the temperature range (0.7 T c — T c ), and reveal a smooth evolution 
across T c . In all cases, the theoretical spectra are calculated for a homogeneous 
system in the normal phase at a temperature T (= 1.1T C ) slightly above T c . In 
panels a-c, the bound (left) and continuum (right) peaks of the experimental 
spectra were independently normalized to the corresponding peak heights, while 
in panel d a single normalization was used . We have adopted the same procedure 
for the comparison. In addition, our calculation provides for the relative weight of 
the bound peak the values 59% (panel a), 44% (panel b), and 33% (panel c). In 
panel a the discrepancy between the positions of the experimental and theoretical 
bound peaks can be accounted for by a molecular multi-channel calculation 22 . 
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Figure 3: Procedure for extracting A m from the tail of the RF spec- 
tra, a, Theoretical RF spectra normalized to the weight A of the continuum 
peak and multiplied by (uj/Ef) 3 ^ 2 - The blue full curve is for (fc^a/) -1 = 2.6 
(with A = 56%) and the red dotted curve for {kpaj)^ 1 = 3.2 (with A = 63%), 
while (kpCLi)^ 1 = and T = 1.1T C in both cases. The horizontal line corre- 
sponds to the value = 0.73 £7? obtained at the same temperature, b, An 
analogous procedure is applied to the corresponding experimental spectra, 
by normalizing the area of the continuum peak to unity. Blue empty squares 
are from Fig. 2b of ref. 7 with (/c^a/) -1 = 2.6 and red full dots are from 
Supplementary Fig. 1 of ref. 7 with (fcpa/) -1 = 3.2, while (kpcii)' 1 = 
and the temperature is below (but close to) T c in both cases. The data yield 
approximately the value 0.25 ±0.10 for the height of the plateau, from which 
we get Aao/Ep = 0.69^0;^ with no input from theory. 
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Figure 4: Diagrammatic representation of the spin correlation func- 
tion above T c . a, DOS contribution where the fermionic f3-f3 single-particle 
propagator (upper line) is dressed by the self-energy E^, which includes pair- 
ing fluctuations of the initial state /3 with its mate a through the pairing 
propagator T a p. b, AL contribution which includes, in addition, pairing 
fluctuations of the final state 7 with the mate a left behind, c, Self-energy 
Eg and pairing propagator Y a p between spins a and (3 in the initial state 
(r a7 between spins a and 7 in the final state is similarly obtained). Solid 
and dashed lines stand for fermionic single-particle propagators and inter- 
actions, respectively, while labels attached to the end points identify the 
relevant spins. All diagrams are drawn in four-momentum space. Only the 
single-particle propagator (5-(5 in panel a is dressed by self-energy insertions. 
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Figure 5: Supplementary Figure 1 Comparison between the DOS 
only and DOS plus AL contributions to the RF spectra on an ab- 
solute scale. The DOS only (red dotted curve) and the DOS plus AL (blue 
full curve) contributions to the RF spectra are compared with each other on 
an absolute scale, defined by the total area of each curve being unity. The 
values of the couplings and temperature correspond to Fig. 2b of the text, 
where the two curves were instead normalized so that the heights of the cor- 
responding continuum peaks coincide with each other. The present absolute 
comparison shows an even more marked difference between the two cases, 
when final-state effects are either neglected (DOS only) or included (DOS 
plus AL). 
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